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Abstract 

We consider a polymer of length N translocating through 
a narrow pore in the absence of external fields. Charac- 
terization of its purportedly anomalous dynamics has so 
far remained incomplete. We show that the polymer dy- 
namics is anomalous until the Rouse time r R ~ N 1+2lJ , 
with a mean square displacement through the pore con- 
sistent with t( 1 +^)/( 1 + 2 ^) ? w ith v « 0.588 the Flory ex- 
ponent. This is shown to be directly related to a decay 
in time of the excess monomer density near the pore as 
£-(i+i/)/(i+2i/) exp(— t/r^). Beyond the Rouse time translo- 
cation becomes diffusive. In consequence of this, the 
dwell-time r^, the time a translocating polymer typically 
spends within the pore, scales as N 2+u , in contrast to 
previous claims. 

PACS Numbers: 36.20.-r, 82.35.Lr, 87.15.Aa 

Transport of molecules across cell membranes is an essential mechanism 
for life processes. These molecules are often long and flexible, and the pores 
in the membranes are too narrow to allow them to pass through as a single 
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unit. In such circumstances, the passage of a molecule through the pore — 
i.e. its translocation — proceeds through a random process in which polymer 
segments sequentially move through the pore. DNA, RNA and proteins are 
naturally occurring long molecules [T] subject to translocation in a variety of 
biological processes. Translocation is used in gene therapy [2], in delivery of 
drug molecules to their activation sites [3], and as an efficient means of single 
molecule sequencing of DNA and RNA [4]. Understandably, the process of 
translocation has been an active topic of current research: both because it is 
an essential ingredient in many biological processes and for its relevance in 
practical applications. 

Translocation is a complicated process in living organisms — its dynamics 
may be strongly influenced by factors like the presence of chaperon molecules, 
pH values, chemical potential gradients, assisting molecular motors etc. [5]. 
In studies of translocation as a biophysical process, the polymer is simplified 
to a sequentially connected string of N monomers. Herein, the quantities 
of interest are the typical time for the polymer to leave a confining cell 
or vesicle, the "escape time" [6], and the typical time the polymer spends 
in the pore or "dwell time" [7], as a function of chain length N and other 
parameters like membrane thickness, membrane adsorption, electrochemical 
potential gradient, etc. [HJ. These have been measured directly in numerous 
experiments |9j. 



Figure 1: Pictorial representation of a translocation event, with the polymer 
shown before, during and after translocation. We number the monomers, 
starting with the end-monomer on the side it moves to. The number of the 
monomer located in the middle of the pore is s. 

Experimentally, the most studied quantity is the dwell time r^, i.e., the 
pore blockade time for a translocation event (see Fig. [I]). For theoretical 
descriptions of r^, during the last decade a number of mean-field type the- 




2 



ories (6j d [8] have been proposed, in which translocation is described by a 
Fokker-Planck equation for first-passage over an entropic barrier in terms 
of a single "reaction coordinate" s. Here s is the number of the monomer 
threaded at the pore (s = 1,...,N). These theories apply under the as- 
sumption that translocation is slower than the equilibration time-scale of 
the entire polymer, which is likely for high pore friction. In Ref. [TO], this 
assumption was questioned, and the authors found that for a self-avoiding 
polymer performing Rouse dynamics, > tr, the Rouse time. Using sim- 
ulation data in 2D, they suggested that the inequality may actually be an 
equality, i.e., ~ tr ~ N 1+2u , which is N 2 ' 5 in two dimensions. Numerical 
data in support of this suggestion in 2D appeared in Ref. [UJ. However, in 
a publication due to two of us, in 3D was numerically found to scale as 
~ N 2AO±om [13], significantly larger than N l+2v ', which is N 2A8 in three di- 
mensions. Additionally, in a recent publication [14] was numerically found 
to scale as i\r 2,52±0 04 in three dimensions (a discussion on the theory of Ref. 
[T4] appears at the end of this paper). Note that these simulations do not 
incorporate hydrodynamical interactions, which are certainly important in 
experiments. Also, these simulations (and all theoretical studies, including 
this one) ignore interactions with other polymers, i.e., they consider poly- 
mers in infinitely dilute solutions, while in cell environments, the solution is 
not infinitely dilute. In this paper we consider translocation in the absence 
of hydrodynamical interactions, and at the end we reflect on the results we 
expect when the hydrodynamical interactions are included. We also note 
here that simulations with hydrodynamical interactions are non-trivial and 
costly. 

Amid all the above results on r d mutually differing by ~ 0(A°- 2 ), the 
only consensus that survives is that > tr [10], [13]. Simulation results alone 
cannot determine the scaling of r^: different groups use different polymer 
models with widely different criteria for convergence for scaling results, and 
as a consequence, settling differences of ~ 0(A 0,2 ) in O(tr), is extremely 
delicate. 

An alternative approach that can potentially settle the issue of Td scaling 
with N is to analyze the dynamics of translocation at a microscopic level. 
Indeed, the lower limit tr for implies that the dynamics of translocation 
is anomalous [10]. We know of only two published studies on the anomalous 
dynamics of translocation, both using a fractional Fokker-Planck equation 
(FFPE) [12], [Hj. However, whether the assumptions underlying a FFPE 
apply for polymer translocation are not clear. Additionally, none of the 
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studies used FFPE for the purpose of determining the scaling of r d . In view 
of the above, such a potential clearly has not been thoroughly exploited. 

The purpose of this paper is to report the characteristics of the anoma- 
lous dynamics of translocation, derived from the microscopic dynamics of the 
polymer, and the scaling of r d obtained therefrom. Translocation proceeds 
via the exchange of monomers through the pore: imagine a situation when a 
monomer from the left of the membrane translocates to the right. This pro- 
cess increases the monomer density in the right neighbourhood of the pore, 
and simultaneously reduces the monomer density in the left neighbourhood 
of the pore. The local enhancement in the monomer density on the right 
of the pore takes a finite time to dissipate away from the membrane along 
the backbone of the polymer (similarly for replenishing monomer density on 
the left neighbourhood of the pore). The imbalance in the monomer den- 
sities between the two local neighbourhoods of the pore during this time 
implies that there is an enhanced chance of the translocated monomer to 
return to the left of the membrane, thereby giving rise to memory effects, 
and consequently, rendering the translocation dynamics subdiffusive. More 
quantitatively, the excess monomer density (or the lack of it) in the vicinity 
of the pore manifests itself in reduced (or increased) chain tension around the 
pore, creating an imbalance of chain tension across the pore (we note here 
that the chain tension at the pore acts as monomeric chemical potential, and 
from now on we use both terms interchangeably). We use well-known laws 
of polymer physics to show that in time the imbalance in the chain tension 
across the pore relaxes as £~( 1 +")/( 1+2i v exp(— £/tr). (Strictly speaking, tr in 
this expression should be replaced by the characteristic equilibration time of 
a tethered polymer with length of 0(N); since both scale as N 1+2u , we use t r 
here, favouring notational simplicity). This results in translocation dynamics 
being subdiffusive for t < tr, with the mean-square displacement (As 2 (t)) 
of the reaction coordinate s(t) increasing as 

£(1+^/(1+2*). diffusiye for 

t>r R . With ^{As 2 {r d )) ~ N, this leads to r d ~ N 2+u . 

We substantiate our theoretical derivations with extensive Monte Carlo 
simulations, in which the polymer performs single-monomer moves. The 
definition of time is such that single-monomer moves along the polymer's 
contour are attempted at a fixed rate of unity, while moves that change the 
polymer's contour are attempted ten times less often. Details of our self- 
avoiding polymer model in 3D can be found in Refs. [HI H6j . 

The key step in quantitatively formulating the anomalous dynamics of 
translocation is the following observation: a translocating polymer com- 
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prises of two polymer segments tethered at opposite ends of the pore that 
are able to exchange monomers between them through the pore; so each 
acts as a reservoir of monomers for the other. The velocity of translocation 
v{t) = s(t), representing monomer current, responds to (pit), the imbalance 
in the monomeric chemical potential across the pore acting as "voltage". Si- 
multaneously, (j)it) also adjusts in response to v it). In the presence of mem- 
ory effects, they are related to each other by 0(t) = j Q dt'fi(t — t')v(t f ) via 
the memory kernel which can be thought of as the (time-dependent) 
'impedance' of the system. Supposing a zero-current equilibrium condition 
at time 0, this relation can be inverted to obtain v(t) = f Q dt'a[t — 
where a(t) can be thought of as the 'admittance'. In the Laplace trans- 
form language, ft{k) = a~ 1 (/c), where k is the Laplace variable representing 
inverse time. Via the fluctuation-dissipation theorem, they are related to 
the respective autocorrelation functions as /i(t — t') = (4>(t)(fi(t')) v=0 and 
a(t - t') = (v{t)vit'))^ =0 . 

The behaviour of n{t) may be obtained by considering the polymer seg- 
ment on one side of the membrane only, say the right, with a sudden introduc- 
tion of p extra monomers at the pore, corresponding to impulse current v (t) = 
p5(t). We then ask for the time-evolution of the mean response (<5$^(t)), 
where <5$^(t) is the shift in chemical potential for the right segment of the 
polymer at the pore. This means that for the translocation problem (with 
both right and left segments), we would have 4>{t) = 6& r '(t) — 5<3>^(t), where 
5& l \t) is the shift in chemical potential for the left segment at the pore due 
to an opposite input current to it. 

We now argue that this mean response, and hence /x(t), takes the form 
/z(i) ~ t~ a exp(— t/r R ). The terminal exponential decay exp(— t/r R ) is ex- 
pected from the relaxation dynamics of the entire right segment of the poly- 
mer with one end tethered at the pore [16J. To understand the physics 
behind the exponent a, we use the well-established result for the relaxation 
time t n for n self-avoiding Rouse monomers scaling as t n ~ n 1+2v . Based 
on the expression of t n , we anticipate that by time t the extra monomers 
will be well equilibrated across the inner part of the chain up to n t ~ 
fX/(i+2v) monomers from the pore, but not significantly further. This in- 
ternally equilibrated section of n t + p monomers extends only r(n t ) ~ n^, 
less than its equilibrated value (n t +p) u , because the larger scale confor- 
mation has yet to adjust: the corresponding compressive force from these 
n t + p monomers is expected by standard polymer scaling [18] to follow 
f/{k B T) ~ 5r{n t )/r 2 {n t ) ~ vp/ [n t r{n t )} ~ t~^ +v ^ +2u \ for p < n t . This 
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Figure 2: Simulation results for the average chain tension component perpen- 
dicular to the membrane proxied by (z^^oo) — z^ 5 >(t)) following monomer 
injection at the pore corresponding to v(t) = p5(t), with p — 10. See text 
for details. Red circles: N/2 = 50, green circles: N/2 = 100, blue circles: 
N/2 = 150, solid black line: r^A 1 ^) w jth v = 0.588 for self-avoiding 
polymers. To obtain a data collapse, the horizontal and vertical axes are 
scaled by (N/2) 1+2u and (N/2) 1+u , respectively. The steeper drop at large 
times correspond to the exponential decay exp(— t/r R ). 



force / must be transmitted to the membrane, through a combination of 
decreased tension at the pore and increased incidence of other membrane 
contacts. The fraction borne by reducing chain tension at the pore leads us 
to the inequality a > (1 + v)/(l + 2u), which is significantly different from 
(but compatible with) the value ot\ = 2/(1 + 2u) required to obtain ~ tr. 
It seems unlikely that the adjustment at the membrane should be dispro- 
portionately distributed between the chain tension at the pore and other 
membrane contacts, leading to the expectation that the inequality above is 
actually an equality. 
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We have confirmed this picture by measuring the impedance response 
through simulations. In Ref. [32], two of us have shown that the centre-of- 
mass of the first few monomers is an excellent proxy for chain tension at the 
pore and we assume here that this further serves as a proxy for 5$. Based on 
this idea, we track (6& r >(t)) by measuring the distance of the average centre- 
of-mass of the first 5 monomers from the membrane, (z^it)), in response 
to the injection of extra monomers near the pore at time 0. Specifically 
we consider the equilibrated right segment of the polymer, of length N/2 — 
10 (with one end tethered at the pore), adding 10 extra monomers at the 
tethered end of the right segment at time 0, corresponding top = 10, bringing 
its length up to N/2. Using the proxy (z^(t)) we then track (5& r \t)). The 
clear agreement between the exponent obtained from the simulation results 
with the theoretical prediction of a = (1 + z/)/(l + 2v) can be seen in Fig. El 
Note that the sharp deviation of the data from the power law £~t 1+l/ V( 1 + 2! -') 
at long times is due to the asymptotic exponential decay as exp(— £/tr), as 
the data collapse shows. 

Having thus shown that fi(t) ~ t~ T +^" exp(— t/r R ), we can expect that 
the translocation dynamics is anomalous for t < Tr, in the sense that the 
mean-square displacement of the monomers through the pore, (As 2 (t)) ~ t 13 
for some (3 < 1 and time t < tr, whilst beyond the Rouse time it becomes 
simply diffusive. The value j3 = a = follows trivially by expressing 
(As 2 (t)) in terms of (translocative) velocity correlations (v(t)v(t')), which 
(by the Fluctuation Dissipation theorem) are given in terms of the time 
dependent admittance a(t — £'), and hence inversely in terms of the corre- 
sponding impedance. 

Indeed, as shown in Fig. [3l a double-logarithmic plot of (As 2 (t)) vs. t is 
consistent with (As 2 (t)) ~ i( 1 +")/( 1 + 2 ^) . The behaviour of (As 2 (t)) at short 
times is an artifact of our model: at short times reptation moves dominate, 
leading to a transport mechanism for "stored lengths" [E] along the polymer's 
contour in which individual units of stored length cannot pass each other. 
As a result, the dynamics of s(t), governed by the movement of stored length 
units across the pore, is equivalent to a process known as "single-file diffusion" 
on a line, characterized by the scaling (As 2 (t)) ~ t 1 ! 2 (not shown here). At 
long times the polymer tails will relax, leading to {As 2 if)) ~ t for t > tr. The 
presence of two crossovers, the first one from (As 2 (t)) ~ t x l 2 to (As 2 (t)) ~ 
t (i+v)/(i+2v) and the second one from (As 2 (t)) ~ £(i+")/(i+2") to (As 2 (t)) ~ t 

at t p& tr, complicates the precise numerical verification of the exponent 
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Figure 3: Double-logarithmic plot of the mean squared displacement of the 
reaction coordinate (As 2 (t)) as a function of time t, for N = 100 (orange), 
200 (red) and 500 (blue). The thick black line indicates the theoretically 
expected slope corresponding to (As 2 (t)) ~ £( 1 +^)/( 1 + 2l/ ). The dashed black 
line corresponds to (As 2 (t)) ~ t 2 ^ 1+2u \ which would have been the slope of 
the (As 2 (t)) vs. t curve in a double-logarithmic plot, if r d were to scale as 
r R ~ N^ 2 \ 

(1 + z/)/(l + 2u). However, as shown in Fig. [3l there is an extended regime 
in time at which the quantity t~^ 1+u ^'^ +2v ' (As 2 (t)) is nearly constant. 

The subdiffusive behaviour {As 2 (t)) ~ t T +^ for t < tr, combined with the 
diffusive behaviour for t > tr leads to the dwell time scaling as ~ N 2+u , 
based on the criterion that \J (As 2 (Td)} ~ N. The dwell time exponent 
2 + v ~ 2.59 is in acceptable agreement with the two numerical results on 
in 3D as mentioned in the introduction of this paper, and in Table I below 
we present new high-precision simulation data in support of r d ~ 111 
terms of the median unthreading time. The unthreading time r u is defined 
as the time for the polymer to leave the pore with s(t = 0) = N/2 and the 
two polymer segments equilibrated at t = 0. Both r u and scale the same 
way, since t u < r d < 2r u [U]. 
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N 




T U /N 2+V 


100 


65136 


0.434 


150 


183423 


0.428 


200 


393245 


0.436 


250 


714619 


0.445 


300 


1133948 


0.440 


400 


2369379 


0.437 


500 


4160669 


0.431 



Table I: Median unthreading time over 1,024 runs for each N. 

We now reflect on the theory presented in Ref. [T4] . 

We have defined r d as the pore-blockade time in experiments; i.e., if we 
define a state of the polymer with s(t) = as '0' (polymer just detached from 
the pore on one side), and with s(t) = N as 'N', then r d is the first passage 
time required to travel from state to state N without possible reoccurances 
of state 0. In Ref. [T4], the authors attach a bead at the s = end of 
the polymer, preventing it from leaving the pore, creating a situation where 
the polymer returns to state multiple number of times before it eventually 
reaches state N. The repeated returns to state implies that by construction 
of the problem, the polymer encounters a free energy barrier on its way from 
state to s = N/2, where the polymer's configurational entropy is the lowest. 
The authors then proceed to express their translocation time (r t hereafter), 
defined as the first passage time required to travel from state to state N 
with reoccurances of state 0, in terms of this free energy barrier. Below we 
settle the differences between r t of Ref. [T4] and our r d . 

Consider the case where we attach a bead at s = and another at s = N, 
preventing it from leaving the pore. We then characterize the state of the 
polymer as follows. At state x and x' the polymer can have all values of s 
except 0,N/2 and N; and at states m and m', s = N/2. The notational 
distinction between primed and unprimed states is that a primed state can 
occur only between two consecutive states 0, or between two consecutive 
states N, while an unprimed state occurs only between state and state N. 
Its dynamics is then given by the sequence of states, e.g., 

, * s 

... N x m x x'O x'm'x'm'x'O x'OxmxmxmxNx'N... 

s V ' 

Td 



9 



where the corresponding times taken (r t and r d ) are indicated. Note in the 
above definitions that r t > T d : since, due to the presence of the entropic 
barrier as described above, r t includes the extra time spent in between the 
first and the last occurrence of state before the polymer eventually proceeds 
to state N. In other words, Tt includes the effect of the entropic barrier, while 
r d does not. A probability argument then leads us to 

n = i = /x(i + /m) (1) 

Td Px + Pm (Pm + Pm')/m(l + fx) ' 

where p m , p m > and p x are the probabilities of the corresponding states, f m = 
Pm/Pm' and / x = p m /p x . The partition sum of a polymer of length n with 
one end tethered on a membrane is given by Z n ~ A n n 71_1 with A a non- 
universal constant and 71 = 0.68 [20], and therefore we have p m + p m > = 

Z 2 n/ 2 / [J2s=o Z sZn- s ] ~ l/N. Similarly, / x ~ l/N [T3]. Finally, f m « 1 [HJ 
yields r t ~ r d . 

In Ref. [U] the authors include a factor A^ 1-71 in r t to account for the 
effect of the entropic barrier. However, we have shown above that r t ~ r^, 
i.e., the free energy barrier does not play a role for the scaling behaviour of 
Tt with N. This implies, since r t includes the effect of the entropic barrier 
and r d does not, that the theoretical expression for T t in Ref. [Hj cannot be 
correct. The numerical result T t ~ _/v 2 ' 52±0 ' 04 in Ref. [14], however, confirms 
our theoretical expression T d ~ N 2+u . 

To conclude, in this paper we have characterized the anomalous dynamics 
of unbiased translocation and obtained the scaling of the dwell time in terms 
of the polymer length. In future work, we will study the role of hydrodynam- 
ics. Rouse friction may be an appropriate model for the dynamics of long 
biopolymers in the environment within living cells, if it is sufficiently gel-like 
to support screened hydrodynamics on the timescale of their configurational 
relaxation. However, we should also ask what is expected in the other ex- 
treme of unscreened (Zimm) hydrodynamics. For our theoretical discussion 
the key difference is that, instead of the Rouse time tr, in the Zimm case the 
configurational relaxation times scale with according to r zimm ~ 

N 3» 

111 

3D, which upon substitution into our earlier argument would gives the lower 
bound value a = (1 + u)/ (3u) for the time exponent of the impedance, lead- 
ing to T d ~ N 1+2u (whose resemblance to the Rouse time is a coincidence - 
note that with hydrodynamics Rouse time loses all relevance). These results, 
however, do need to be verified by simulations incorporating hydrodynamics. 
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